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fH Abstract 

In a system made up of independent random walks, fluctuations of order n^/^ from 
the hydrodynamic limit come from particle current across characteristics. We show 
i_i that a two-paxameter space-time particle current process converges to a two-parameter 

Gaussian process. These Gaussian processes also appear as the limit for the one- 
dimensional random average process. The flnal section of this paper looks at large 
deviations of the current process. 
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It is well known that particle systems that appear different at the microscopic level often 
behave almost identically at a macroscopic level. This has been observed in the hydrody- 
namic limits and fluctuation results of several particle models. Consequently, there is much 
to be gained in studying the behavior of simpler stochastic particle systems in the hope 
that at the macroscopic level they will reflect the behavior of a universal class of systems. 
While the hydrodynamic limit of several models have been studied, fluctuation results have 
proved elusive for many systems. In this paper we consider particle current fluctuations in 
the one dimensional independent random walk model. 
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The hydrodynamic limit for particle distribution in typical asymmetric systems are 
solutions to p.d.e's of the form 

dtu + d,f{u) = o. (1.1) 

In the case of nearest neighbor Totally Asymmetric Simple Exclusion Process (TASEP) 
in one dimension, the flux function /(p) = p(l — p). For non-interacting particle systems 
/(/?) = V ■ p where v is the average velocity of the particles. Thus the relevant p.d.e. for a 
system of independent asymmetric random walks is 



dtu + vd^u = (1.2) 



(Prop 3.1, page 15 in [5]) 



From the transport equation (1.2) we see that in the independent random walk model. 



the initial density profile shifts with velocity v. Consider an observer moving at constant 



velocity v. The path of the observer is a characteristic line of the transport equation (1.2). 
It is natural to expect the net current of particles across the path of the observer to be zero. 
But what are the fluctuations in this particle current? This is the question we address in 
this paper. It has been observed that these current fluctuations are of order n^^^ [13]. Here, 
n is the scaling parameter. Typically we scale both space and time by n in asymmetric 
models, this is called Euler scaling. In symmetric models we use diffusive scaling i.e. we 
scale time by n and space by ^/n. There is a general belief that when / is linear (i.e. 



/" = 0) in (1.1), the fluctuations in particle current across characteristics of (1.1) should 
be of order n^. This has been shown for the random average process (RAP) and for the 
one dimensional independent random walk model where /" = 0. 

In this paper we study both the fluctuations and the large deviations of the current 
process for independent walks. For the fluctuations we consider the current process indexed 
both by time and spatial shifts of order ^/n of characteristic lines. The y/n order for spatial 
scaling is the natural one because the individual random walks fluctuate on that scale. We 
extend the distributional limit of [13] to a process limit for the space-time current process. 
The space-time current process was also studied for RAP in |T] but only convergence of 
finite dimensional distributions was shown without process-level tightness. The same family 
of Gaussian processes arises as limits for both RAP and independent random walks. 

It is interesting to note that there are models which are not asymmetric yet exhibit 
subdiffusive current fluctuations with Gaussian scaling limits. It was conjectured (conjec- 
ture 6.5 in [l3]) that subdiffusive fluctuations in 1 dimensional nearest neighbor symmetric 
simple exclusion processes (SSEP) converge to IBM with Hurst parameter 1/4. It was sub- 
sequently proved in the finite dimensional distributions sense and has recently been proved 
in the full functional central limit theorem sense in [11'. This says that the universality class 
of current fluctuations of order n^/^ contains both symmetric and asymmetric processes. 
However, the symmetric and asymmetric processes differ on the level of hydrodynamics. 
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This paper is organized as follows. We start with a description of the independent 
random walk model and the statements of the main results in section |2j The next three 
sections |3] |4] and [5] cover the proofs. Section [3] gives the convergence of finite dimensional 
distributions and section |4] proves process level tightness. A note on the tightness methods 
used here: since we are interested in a two-parameter process, the standard theorems on 
convergence in L']r[0,oo) and C]r[0,oo) spaces do not apply. We appeal to two papers, [2] 
and [6], that provide suitable criteria for deducing tightness. [2] gives the context in which 
we speak of convergence for two-parameter processes and a tightness criterion. The proof 
of Proposition 5.7 in [6] is extended to two dimensions to prove the tightness criterion. The 
last section contains proofs of some large deviation results for the current process. 

2 Model and results 

2.1 Independent random walk model 

Consider particles distributed over the one dimensional integer lattice which evolve like 
independent continuous-time random walks. We are given the initial occupation variables 
Vo = ivoix) : X £ Z} defined on some probability space. Let Xm,j{t) denote the position at 
time t of the jth random walk starting at site m. The common jump rates of the random 
walks are given by a probability kernel {p{x) : x G Z}. Once the initial positions of the 
random walks are specified, their subsequent evolutions {Xm,j{t) — Xm,j{0) : m £ j = 
1, . . . , r/o(?Ti)} are as i.i.d. random walks starting at the origin, on the same probability 
space, independent of rjQ. Define 

meZ j=l 

to be number of particles on site x at time t. 

Assumption 2.1. For the random walk kernel, we assume that, for some 5 > 

^e'^^'pix) <oo for \e\< 5 (2.1) 

(This assumption will enable us to calculate large deviation hounds for the random walks.) 

Throughout this paper, we assume that N denotes the set of positive integers. Let 
{r/Q : n e N} be a sequence of initial occupation variables defined on some probability 
space. 
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Assumption 2.2. For each n, the initial occupation variables {r?o(^) '■ x ^ l^} are inde- 
pendent. They have a uniformly hounded twelfth moment: 

sup E[n'^{xY'^\ < oo. (2.2) 

Let Pq{x) = Etjq^x) and Vq{x) = Var[r]Q{x)] be the mean and variance resp. of the initial 
occupation variable r/Q(x),x G Z. Let po and vq be two given nonnegative, finite numbers. 
The means Pq and variances Vq approximate po and vq in the following sense: There exist 
positive integers L = L{n) such that n~^/^L{n) and for any finite constant A, 



lim sup n^^^ 

\m\<Ay'n log n 



L(n) 



-^^pS(m + j)-Po|=0 (2.3) 
^"^f j=i 



The same assumption holds when p^ and po are replaced by Vq and vq. 



As in |I3], the reason for tlie complicated assumption (2.3 1 is to accommodate both 
random and deterministic initial conditions. For random rj^^x) we could take Pq{x) = 
po(^) for some sufficiently regular function po{-). However, for deterministic rjQ^x) we 
cannot do this unless po{x) is integer- valued. A couple of examples illustrating random 



and deterministic initial configurations that satisfy assumptions |2.l| and 2.2 can be found 
in [13]. 
Let 

V = j;p(x) and K2 = ^^x^p(x). 

X X 

The characteristics of ( |1.2[ ) are straight lines with slope v. Fix T > and S > 0. For 
t £ [0, T] and r G [—5, S], we let Yn{t, r) denote the net right-to-left particle current during 
time [0, nt] across the characteristic line starting at ([r-y/n],0). More precisely, 



oo Voim) 

Yn{t,r) = J2 Yl [HXm,j{nt) < [nvt] + [rV^]}l{m > [r^^]} 

m=— oo j=l 

- l{Xmjint) > [nvt] + [r^]}l{m < [r^n]}] 



(2.4) 



where Xkj{-) is the jth random walk that starts at site k. Note that the random walks 
denoted as Xkj in the definition of Yn{t,r) should actually be X]^^, but we drop the 
superscript n for notational simplicity. 

2.2 Distributional limit 

We give a brief description of the path space of the process Yn{-,-). Let D2 = D2{[0,T] x 
[— 5, 5],M) be the space of 2-parameter cadlag functions with Skorohod's topology. Let 
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Q := [0,T] X [—5,5']. For any (t,r) G Q, we can divide Q into four quadrants: 

Qlt,r) ■= q) (^Q ■■ s>t,q>r}, Qj^ := {(s, G Q : s > t, g < r}, 

Qft,r) e Q : s < g < r}, Qj^ := {(s, G Q : s < t, g > r}. 

Then the precise definition of D2 is 

D2 = {f : Q ^ ^ such that for any point (t, r) G Q, hm /(s, q) 

exists for i = 1, 2, 3, 4 and hm f{s,q) = f{t,r)}. 

{^,g)6Qjt,,) 

(s,g)^(t,r) 

In other words, D2 contains functions that are continuous from the right and above with 
hmits from the left and below. Skorohod's topology in Z?m[0, 00) is extended to this space in 
the most natural way. The space of multiparameter cadlag functions and their topology is 
described in detail in [2] . By Theorem 2 in [2] , a sufficient criterion for the weak convergence 
Xn ^ X in D2 is, 

1. For all finite subsets {{ti,ri)} C [0,r] x [-S,S], 

{Xniti,ri), - ■ ■ ,Xn{tN,rN)) (X(ti,ri),--- ,X{tN,rN)) 

weakly, and 

2. lim^^o li™sup„ > e} = for all e > 0, where the modulus of continuity is 
defined by 

Wx{6) = sup \x{s,q) - x{t,r)\. 

is,q},(t,r)e[0,T]xl~S,S] 
\is,q)-(t,r)\<5 

Clearly, {Yn{t, r) : t £ [0, T], r G [—5, S]} are D2-valued processes and we can use the above 
criterion to prove their convergence in the weak sense. 

Denote the centered Gaussian density and distribution with variance cj^ by 

11 

Define also 

^^2{x) = a^cp^x) - x{l - $^2(x)). 
For {s,q),{t,r) G [0, T] x [—5", 5], define two covariances by 

Toiis, q), {t, r)) = ^,,s{\q - r\) + ^.,t{\q - r\) - - r|) (2.5) 

and 

Tg{{s,q), (t,r)) = ^^,(^t+s){\q " r\) - ^.,it-s\{\q " r\). (2.6) 
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THEOREM 2.1. Define Yn{t,r) as in (2^. Then under Assumptions \2l\ and \K2 



as n ^ oo, the process n •) converges weakly on the space D2 to the mean zero 

Gaussian process Z{-,-) with covariance 

EZ{s, q)Z{t, r) = poTg{{s, q), {t, r)) + VoTo{{s, q), {t, r)). (2.7) 

Note: We will show later that n-^/^EYn{t,r) as n ^ 00 uniformly for t € [0, T] 
and r G [— S, S"]. Hence, in the above theorem we can replace n~^^'^Yn with the centered 
current process with impunity. 

COROLLARY 2.1. Under the invariant distribution where {^"(x) : x G Z} are i.i.d. 

Poisson with mean p for all n, n~^^^Yn{-, •) converges weakly in D2 to a mean-zero Gaussian 
process Z{-,-) with covariance 

EZ{s,q)Z{t,r) = p{^.U\q " ^1) + "^^AW " r\) - ^.,|t-.|(k " H))- 
In particular, for a fixed r the process {Z{t,r) : t G [0,T]} has covariance 



EZ{s, r)Z{t, r)=p^^{^s + Vi- V\t^) , 
i.e., process Z(-,r) is fractional Brownian motion with Hurst parameter 1/4. 



The covariance (2.7) is the same as the covariance of the limiting Gaussian process in 
the RAP model found in |T], with different coefficients in front of Tq and Tq. As in [Ij, we 
can represent the Gaussian process Z{-,-) as a stochastic integral: 

Z{t, r) = y/K2Po / 4>K2(t-s) {r - z)dW{s, z) 

+ sgn{x - r)^i^^t{-\x - r\)dB{x). 



The equality in (2.8) is equality in distribution. is a two-parameter Brownian motion 
defined on M_|_ x M, i3 is a one-parameter Brownian motion defined on M, and W and B 
are independent of each other. The stochastic integral clearly delineates the two sources of 
fluctuations in the current. The first integral represents the space-time noise created by the 
dynamics, and the second integral represents the initial noise propagated by the evolution. 

2.3 Large deviation results 

We first state large deviation results for Yn{t, r) with fixed r G M and t > 0. 
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Assumption 2.3. Assume the initial occupation variables {'r]Q{m) : m G Z} are i.i.d. Let 

7(^) = logEe'^''oM 

with effective domain 



-D-y := {a G M : 7(a) < 00}. 



Assume Dj = M. 
For A G M, define 

and 

with effective domain 



log{$«2t(y) + e^(l - for y > 

log{e-^$«2t(y) + 1 - for y < 



/oo 
7{Zx{y))dy 
-00 



-Da := {a e M : A(a) < 00}. 



A turns out to be the Hmiting moment generating function of the current. 
Recall 

Definition 2.1. A convex function A : M — ^ (—00,00] is essentially smooth if: 

a) D"^(interior of D\) is non-empty. 

b) A(-) is differentiable throughout D\. 

c) A(-) is steep, i.e., lim„_^.oo |VA(A„)| = 00 whenever {A^} is a sequence in converging 
to a boundary point of D\. 

Let 7*(x) := sup;^g]g{A • x — 7(A)} be the convex dual of 7(-). Let 

I{x) := sup{A • X - A(A)}. 



Recall the usual definition of Large Deviation Principle (LDP). The sequence of random 
variables {Xn} satisfies the LDP with rate function J{x) and normalization {\/n} if the 
following are satisfied: 

1. For any closed set F C M, 

lim sup ^= log P{Xn G -F} < — inf J{x) 
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2. For any open set G C 



liminf logP{X„ G G} > - inf J{x) 



The rate function J is said to be a good rate function if its level sets are compact. 

Let Brp{X) := logEe'^^ , where X ~Bernoulli(p), be the logarithmic moment generating 
function of Bernoulli random variables. Its convex dual is 

X 1 — X 

Br*(x) = xlog — h (1 — x) log 

^ p 1 — P 



for < a; < 1. Define a{x) implicitly by 

X = A'(a(x)). 



(2. 



This definition is well defined for all x G M since A(-) is strictly convex. For any a G 
define 

^"(y) ■= 5^ , — for y G 



1 - $«2t(y) + e-°$«2i(2/) 



By (2.9) we have 



Define 



and 



h{x) 



7* {l'{Za(^^){y))} dy 



h{x) := / 7'(Z„(^.)(y))Sr; (j^)(F^(^.)(y))dy. 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



Recall that we assumed (') ^^'^ i.i.d. at the beginning of this section. Consequently, the 
underlying distribution is shift invariant. The rate function in Theorem 2.2 therefore does 
not depend on r, as the marginals of the current process, Yn{t, r), are shift invariant. 



THEOREM 2.2. Let Assumptions 



2.1 



and 



2.3 hold. For fixed real r andt > 0, n"^/^y„(t, r) 



satisfies the large deviation principle with normalization {\/n} and good, strictly convex rate 
function 



I{x) = h{x) + hix), xG 



(2.14) 
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A few words on the rate function. / has a unique zero at zero. The rate function / 
balances two costs: the cost of deviations in the initial occupation variables, given by Ji, 
and the cost of deviations in the probability with which particles cross the characteristic 
lines, given by l2- In the macroscopic picture, 1 — ^K2t{y) is the a priori probability with 
which particles initially at distance y > from the characteristic line cross it by time t, 
while a particle at distance y < crosses the characteristic line with probability ^f^^tiv)- 

An intuitive understanding of the LDP is as follows. To allow a current of size x at 
time t the system deviates in such a way that its behavior is governed by a new probability 
measure. Under this new probability measure the mean number of particles initially at 
site y is 7'(ZQ,(^)(y)) and the probabilities 1 — ^K2t{y) and ^^^tiy) are tilted to give new 
probabilities 1 — F^(^^^{y) and F^(^^-^{y). The term Br^^^^(^y^{F^(^^^{y)) measures the cost of 
the deviation of the probability 1 - ^>K2t(y) (or ^K2t{y)) to 1 - Fa(^){y) (or (?/)). The 
new measure depends on a{x) which is chosen so that the mean current under the new 



measure is x, this is evident from (2.11) 



The modus operandi for the proof involves using non-rigorous, intuitive ideas for finding 
a candidate for the rate function and then checking if it is, in fact, the correct rate function. 

We give explicit formulas for the two simplest cases: the stationary situation with i.i.d. 
Poisson occupations, and the case of deterministic initial occupations. 



COROLLARY 2.2. 1. 

the rate function is 



U Voi') ~ Poisson{p), then under assumptions 2.1 and 2.3 



I{x) = xlog 



f X^pTX 



+ Wi + 




for X G M. 
(2.15) 



2. IfrjQ^-) = 1, then under assumption 2.1 the rate function is 

where Fn(^.\(y) = - — . , ^ , ^l^/ifi — rr with a(x) chosen so that 



dy 



^1 - Fa(x){y))dy 



Fa(x){y)dy. 



For the process level, under the stationary distribution, we can show an abstract LDP 
by applying a theorem from [7] . But currently we do not have an attractive representation 



of the rate function. The rate function is given in terms of a variational expression in ( 5.14 ) 
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THEOREM 2.3. //?7g (•) ~ Poisson{p), then under assumptions 2.1 and 2.3 the sequence 
of processes {n~^/^y„(-, 0)} in L'k[0,oo) satisfies a LDP with a good rate function. 

There exist some large deviation results for independent random walk systems in the 
literature These papers essentially deal with large deviations of occupation times 

of sites. Lee [9] and Cox and Durrett [3 find LDP's for weighted occupation times of 
sites under deterministic and stationary (i.i.d. Poisson) initial distributions. Even the 
normalizations of the LDP's were distinct for these two cases. This is in contrast to our 



current large deviations in Theorem 2.2 where the normalizations for the Poisson and the 
deterministic case were the same. 

3 Weak Convergence of Finite Dimensional Distributions 



We begin the proof of Theorem |2.1| by first showing weak convergence of the finite dimen- 
sional distributions of n^^/^Yn{-, •). 



PROPOSITION 3.1. Define Yn{t,r) as above. Then under assumptions 2.1 and 2.2 the 

finite dimensional distributions of the processes {l^(t,r) : (t, r) G [0,T] x [— S, S*]} converge 
weakly as n ^ oo to the finite- dimensional distributions of the mean zero Gaussian process 



{Z{t,r) : {t,r) G [0, T] x [—5", 5]} with covariance given in (2.1). 



To prove convergence of finite dimensional distributions we use Lindeberg-Feller and 
check the conditions of Lindeberg-Feller by brute force. We show that the expected value 
of the process converges uniformly to and hence we can consider the centered process 
when proving convergence. We also appeal to large deviations of random walks to control 
the contributions to the current process from distant particles. 

Fix N space-time points: (ti, ri), (t2, 7'2)i • • • i {tN-, ^n) where (tj, rj) / {tj^fj) for i 7^ j- 
We will prove that as n — > oo the vector 

n-i/^(y„(ti, ri), y„(t2, ra), ...,Yn{tN, r^)) 

converges in distribution to the mean-zero Gaussian random vector 

(Z(ti,ri), Z(t2,r2), Z{tN,rN)) 

with covariance 

EZ{ti,ri)Z{tj,rj) = poTg{{ti,ri),{tj,rj)) + voro{{ti,ri), {tj,rj)). 
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Let = {9i, On) € R be arbitrary. Recall that Xkj{-) is the jth. random walk that 
starts at site k. 

N 
i=l 

N f CO VoiiTi) 

= ^'^^^ Yl H {l{Xm,j{nti) - Xm,j{0) < [nvti] + [n^M - m}l{m > [n^M} 

i=l ^m=— oo j=l 

- l{Xm,jinti) - Xm,jiO) > [nvti] + [riVn] - m}l{m < [viVn]}) | 
Let 

^mj = {^rnA^t) - X^j{0) < [nvt] + [r^^ - m} 

and 

A^'"^ = {X{nt) < [nvt] + [r^n] - m} 

where X{-) represents a random walk with rates p{x) starting at the origin. The evolution 
of the random walks is independent of their initial occupation numbers ?7o(a;). Define 

C/^(i,r) :=n-V4 ^ {l{<^.}l{m > [ry^} - l{{AX,r}l{m < [ry^]}} 

- n-^/^p^{m){P{A*'''}l{m > [r^M} - P{{A^'''y}Hm < [r^M}} 

and 

N 

Um ■= ^OiUmithn). 
1=1 

Then we can write 

N N N 

1=1 1=1 1=1 

oo N 

= J2 Um + n-^'^Y^iEYn{ti,n). 

m=—oo 1=1 

We split J2m=-oo ^"^0 sums 5i and 52 as follows. Choose r{n) so that r(n) = 

o{^/log n) and r(n) —>■ oo slowly enough that 

r{n)H{n^^^) — as n — oo, 
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where H{M) = sup^>^^^^j^E[r]^{x)H{ri^{x) > M}]. Write 

AT 

n-V4^^,{y„(ti,r,)-£;y„(ti,n)}= ^ Um+ Yl 

i=l \m\<r{n)^/n \m\>r{n)^/n 

=: Si + 52. 

We show that 5*2 goes to in and use Lindeberg-Feller for Si to show that it converges 
to a mean-zero normal distribution. 

LEMMA 3.1. ES^ as n ^ oo. 

Proof. EUjn = 0, so EUm = 0. S2 is therefore a sum of independent mean zero terms Um- 
By Schwarz Inequality, 

N 

ESl= Eul<\\efY E Euliu.n). 

\m\>[r{n)^/n\+l «=1 \m\>[r{n)^]+l 

Since N is fixed, it suffices to show for fixed (i, r) 

00 

lim EUl{t,r) = 0. 

n— >oo ^— ' 

m>[r(n)\/n]+l 

Recall that = {Xm,j{nt) - Xm,j{0) < [nvt] + [r^/n\ - m} and A*'"" = {X{nt) < 

[nvt] + [r^/n] — m}. If Tk = X^^Li Zi is a random sum of i.i.d. summands Zi independent 
of random K, then 

Var[TK] = EK ■ Var Zi + {EZif ■ Var K. 
For fixed m, take K = rfQ{m) and 

Zj = l{<,.}l{m > [rv^} - l{{AXjr}l{m < [rV^]}. 

Then, 

E[U^{t,r)] = Var[n-V4 ^ {l{Al;:jl{m > [rV^]} - l{{A'^^^r}l{m < [rV^}}]. 

Var Z, =P{A'^^^)P{{A'^^^r)l{m > [r^n]} + P{A'^^^)P{{A'^^^r)l{m < [r^n]} 
= P{X{nt) < [nvt] + [r^/n] - m)P{X{nt) > [nvt] + [r^/n] — m). 
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[EZif = P{X{nt) < [nvt] + [ry^] — m) \{m > [r^/n]} 

+ P{X{nt) > [nvt] + [ry/n] - mfl{m < [r^/n\}. 

E[Ul{t,r)] = n-^/^pUm)P{A''')P{{A''Y) + n-^/\i;{m)[P{A''^fl{m > [r^]} 



Using the uniform bound (2.2) on moments we get 



E[Ul,{t,r)] < n-^/^C[P{X{nt) < [nvt] + [r^/n] - m)l{m > [r^/n]} 
+ P{X{nt) > [nvt] + [r^/n] - m)l{m < [r\/n]}] . 

By standard large deviation theory, for arbitrarily small 6, there is a constant < K < 
oo such that when m > [ry^, 



P(X(nt) <[nvt] + [r^/^]-m) < < ; ^, r l v Ji - .g^-^ 



and when ni < [vy/n], 

( g{-K{[r^-m)^/nt} ■£ I U - ml < nt6 

P(X(nt)>M] + [rV^]-m)< \- ,^ (3.2) 



if |[rA/n] — m\ > nt5. 



Consequently, 



E E[Ul{t,r)] 

\m.\ = [r{n)^i]+l 

[ntS\ oo 

mi = [r(n)v'n] — [ry^+l mi = [nt(5]+l 

[nt(5] oo 

mi = [r {n)y/n\ + [r y/n\ +1 mi = [nt(5] + 1 

^OO ^ /"OO ^ 1 

^ _^-KntS ^ 



f'OO f'OO 

<C e-^'^'I'dx + C e-^'^'/'dx + 2 

J r(n)—r J r(n)+r 



lr(n)—r J r{n)+r ^ 

as n — > cxD (since r(n) co). Therefore, -E^l ^ as n ^ oo. □ 
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Define 



and 



/■oo 

rQ>{{s,q),{t,r)) = / P{,/]^B{s) <q-x)P{,/l^B{t) <r-x)dx 

JqVr 

/qAr 
P{^B{s) > q- x)P{y/J^B{t) >r- x)dx 
-oo 

- l{r >q} I P{^B{s) <q- x)P{^B{t) >r- x)dx 

Jq 

- l{q > r} P{^/k^B{s) > q - x)P{^/K^B{t) < r - x)dx 



/oo 
{P{^2B{s) < q - x)P{^2B{t) >r-x) 
-oo 

- Pi^/K^B{s) <q-X, y/K^Bit) > T - x)}dx] 



B{-) is standard Brownian motion. Tlie covariance terms (2.5) and (2.6) appear in the 
calculations as Fq^^ and Fg^^ resp. 

LEMMA 3.2. 5*1 converges to the mean-zero normal distribution with variance 



Proof. We apply Lindeberg- Feller to Si. We check: 
1. For any e > 



lim V E[Ull{\Um\>e}]=G 

n. — >rya ^— • 



(3.3) 



|m|<[r(n)v^ 



and, 



2. 



lim ^[^™] = 

n— >oo ^— ' 

|m|<[r(n)v^] 



(3.4) 



\Um{t,s)\ < n ^^^rjQ^m) + n ^^^pQ^m). By uniform bound on moments (2.2), we get 

\Um\<Cn-'^^[r]^{m) + l]. 
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By choice of r(n) we get (3.3): 



lim V E[Ull{\Um\>e}]=d. 

71. — >no ^ — * 



|m|<[r(n)v^ 

We prove ( |3.4[ ) from the lemma below. 

|m,|<r(n)0i ^<i,j<N \m\<r{n)^ 



Let 



5= E[U„,{s,q)Urn{t,r)]. 

|m|<[r(n)v^ 



LEMMA 3.3. 



lim S = /9or«((s, g), (t, r)) + t^oLf ((s, g), (t, r)). 



Proof. 



E[U^{s,q)Um{t,r)] 

= n-V2cov[x; {l{yl^^ll{m>[gV^}-l{(A^.ni{m<[gV^]}} , 

{l{<,}l{m > [rV^]} - l{{AX,r}l{m < [r^}} 

Assume Zi are i.i.d. random variables independent of the random nonnegative integer 
K, then 



K 



K 



Cov[^ f{zi),Y, 9{zi)] = EK ■ Cov[/(zi), + Var K • Ef{z{) ■ Eg{zi). 



1=1 



i=l 



For fixed m, take K = rjQ^m), 

f{Xm,) = l{A'^',}l{m > [qV^} - liiA^^'^mm < [q^} 

and 
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We now adopt the simpler notation 

A = {X{ns) < [nvs] + [q\/n\ — m} 

and 

B = {X{nt) < [nvt] + [ry/n] - m}. 
A straightforward computation gives 

Cov[/(X™,i),5(X^,i)] = P{A)P{B^) - P{AnB^). 

Next we compute 



Ef{Xm,i)Eg{Xm,i) =P{A)P{B)l{m > [gV^}l{m > [r^^} 

+ P{A'')P{B'')l{m < [q^Jn]}l{m < [r^Jn} 

- P{A)P{B'')l{m > [q^Jn]}l{m < [r^Jn} 

- P{A'')P{B)l{m < [q^Jn]}l{m > [r^Jn}. 

Putting these computations together we get 
S =n-i/2 Yl [pUm){PiA)P{B'^) - P{A n B^)] 

\m\<[r{n)yW\ 

+ vl{m) P{A)P{B)l{m > [qy/n\ V [r^/n\] + P(A'=)P(B^)l{m < [qy/n\ A [ry/n] 

- P{A)P{B'')l{[qy/^ <m< [rVn} - P{A'')P{B)l{[r^/^ <m< [qV^]} 

X{nt) is a sum of Poisson(nt) number of independent jumps .^j, each jump distributed 
according to p{x). Therefore, 

Var(X(nt)) = nt ■ Var(^i) + {E^if ■ (nt) = nt{K2 - v^) + v^nt = ntK2. 

By Donsker's Invariance Principle the process {{X{nt) — [nvt])/y/nK2 '■ t > 0} converges 
weakly to standard 1-dimensional Brownian motion. Therefore, as n — > oo, assumption 



(2.3) and a Riemann sum argument together with the large deviation bounds on X{nt) 
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gives us 

/oo 
{P{B{k2s) <q- x)P{B{K2t) >r-x) 
-oo 

- P{B{k2s) <q-x, B{K2t) > r - x)]dx 

+ Vo\ [ P{B{k2s) <q- x)P{B{K2t) <r- x)dx 

VJqWr 

/qAr 
P{B{k2s) > q- x)P{B{K2t) >r- x)dx 
-oo 

- I [l{q < r}P{B{n2s) <q- x)P{B{K2t) > r - x) 

JqAr 

+ l{r < q}P{B{K2s) > q- x)P{B{K2t) < r - x)] | 

= por«((s, q), {t, r)) + voT^o\is^ Q)^ 0)- 

The reasoning behind the convergence of S to the above hmit is the same as in Lemmas 4.3 
and 4.4 of [13]. The reader is referred to [13 (page 778) for a more detailed explanation. 

□ 



This verifies ( |3.4[ ): 

Jim Yl ^[^-]= E 0,e,{poTW{{U,r,),{t„rj))+voT^;^\iU,r,),{t,,r,))}. 

\m\<r{n)y^ l<i,j<N 

We can now conclude from Lindeberg-Feller that Si converges to mean-zero normal 
distribution with 

a^= Yl OiOj{po^'i\iU,ri),{tj,rj))+voT^o\{U,ri),{tj,rj))}. 

□ 

LEMMA 3.4. 

lim sup n~'^/'^\EYn{t,r)\ =0. 

n^oo 0<t<Tfl<\r\<S 

Proof. The proof is similar to the proof of Lemma 4.6 (page 781-783) in jl3j . □ 
Showing that 

r«((s, q), (t, r)) = ^^,^t+s){\q - r\) - ^ n,i\s-t\){\q - r\) = r,((s, q), (t, r)) 
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and 



is an exercise in calculus. 



T^^\{s,q),{t,r))=To{{s,q),{t,r)) 



Proof of Proposition 3.1[ By Lemmas 



3.1 



3.2 



and 



3.4 



we have, as n — > oo, n J2i^=i ^i^niUi 



converges to a mean-zero normal distribution with variance o"^, where 

^<i,j<N 

for arbitrary (e'l, 6'Ar) G R^. We can therefore conclude that as n ^ oo the vector 
n~^/'^{Yn{ti,ri),Yn{t2, ^2), ...,Yn{tN, i^n)) converges in distribution to the mean-zero Gaus- 
sian random vector {Z{ti,ri), Z{t2,r2), Z(ti\i,ri^)) with covariance 

EZ{ti,ri)Z{tj,rj) = poTq{{ti,ri),{tj,rj)) + voro{{ti,ri), {tj,rj)). 



□ 



4 Tightness and completion of proof of Theorem 2.1 

In this section we first develop a criterion for tightness for processes in D2. The tightness 
criterion is in terms of a modulus of continuity. We then proceed to check if our scaled 
current process satisfies the tightness criterion. The following is an extension of Proposition 
5.7 in [6j to two-parameter processes. WLOG for simplicity we replace the region [0,T] x 
[-S, S] with the unit square [0, 1]^. For any h £ D2 = D2{[0, 1]^M), define 

WhiS) = sup \h{s,q) - h{t,r)\. 

s,t,q,r&[0,l] 
\{s,q)-{t,r)\<5 

PROPOSITION 4.1. Suppose {Xn} is a sequence of random elements of D2 = D2{[0, 1]^,M) 

satisfying these conditions: 

For all n there exists 6n > such that 

1. there exist P > 0, a > 2 , and C > such that for all n sufficiently large 

E{\Xnis,q)-Xnit,r)f) < C\is,q) - (t,r)r (4.1) 
for all s,t,q,r £ [0, 1] with \ (s, q) — (t, r) \ > 5n, and 

2. for every e > and r/ > there exists an no such that 

P{wxA^n) > e) < ri for all n > no. (4.2) 
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Then, for each e > and r] > 0, there exists a 5, < 6 < 1, and an integer uq, such that 

P{wxS^) > e) < ??, for n > uq. 

To prove this proposition we require the following lemma. 

LEMMA 4.1. Let < ko < k. If points {s,q) and {t,r) lie on the grid i.e. s = 
a — j- t—— r—^ nnd\^ — 1<2~'^" I J ^l<2^^o then 

1. it is possible to move from (s, q) to {t, r) in steps of size , ko < h < k, moving one 
co-ordinate at a time, where a step of size 2~^, for any h, occurs at most 4 times; 

2. also, we can choose our steps in such a way that we make a jump of size 2^^ only if 
we lie in the 2~^ grid. 

Proof. We fix /cq and prove the lemma by induction on k. When k = k^, {s,q), {t,r) 
£ 2-''m X 2-'="N. We are given that \A- - ^\ < 2"*^", li - Xrl < 2-^o. So, either 

o I 2 I — 'I 2 2 I — ' 

both points coincide or there is a difference of 2~^° in one co-ordinate or both co-ordinates 
between these points. We can therefore move from (s, q) to {t, r) in at most 2 steps, each 
of size 2"'^", moving one co-ordinate at a time. Clearly condition (2) also holds as the only 
jumps possible here are of size 2"'^" and we lie in the 2~^° grid. 

Let k > ko, {s,q) = (^,^), {t,r) = {^,^) ,| ^ - < 2-'=o and |^ - 1^1 < 2-^<\ 
Either the points (s, q) and (t, r) already lie on the 2~^^~^^ grid, or if not, the points (s, q) 
and {t,r) are each at most two jumps of size 2~^ away from the 2~(^~-^)-grid. We can 
therefore move both {s,q) and (t, r) closer to each other and onto the 2~(^^^''-grid in at 
most four jumps of size 2~^. 

By the induction hypothesis and by our choice of jumps, we are done. □ 

Proof of Proposition \4.1\ We may assume without loss of generality that (5„ = 2"'^ for 
some k = k{n) > depending on n. This is because if we have P{wx„{^n) > e) < rj 
for some (5„ > 0, then we can find a k{n) > with 2^'^^'^^^^ < 5n < 2~^^^^ such that 
«;x„(2"'=(")) < 2wx„{Sn). So P(w;x„(2-'=(«)) > 2e) < P{wxM > e) < r,. Therefore it is 
sufficient to prove the theorem for 6n = 2"*"'*^") for k{n) > 0. 

Given n and 6, ii 6 < 6n, then wx^i^) < wx^i^n)- If < ^, then for any two points 
(s, q), {t, r) with < |t — s|, |g — r| < (5, we have 

\h{s,q) - h{t,r)\ < \h{s,q) - h{s',q')\ + |/i(t,r) - h{t',r')\ + \h{s',q') - /i(t',r')|, 

where s',t',q',r' £ (5„N, < \s' -t'\, \q'-r'\ < 5 and < \t-t'\,\s- s'\, \r-r'\,\q-q'\ < 6n. 
Thus for any 6n < S, 

Wh{5) < 2wh{Sn) + sup \h{s,q) - h{t,r)\. 

t,s,g,rG(5„N)n[0,l] 
|t-s|,|r-ij|<(5 
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Therefore by (4.2), we only need to show that for any e > and rj > 0, there exists a 5 > 

(4.3) 



such that for all n sufficiently large 



where 



wl\d) = sup \h{s,q) - h{t,r)\. 

t,s,g,re(<5„N)n[0,l] 
\t-s\,\r-q\<5 



This follows from (4.1) by the following "dyadic argument": 



Choose a A such that 2^^ '^^ < < 1. Given n for which (4.1) is satisfied, let 

Gk = {\Xni^, ^) - ^n(^, < A'^ for i = 0, 1, . . . , 2'= - 1, 
i = 0,l,...,2^} 

and 

Hk = i) - ^«(^' ^)| < A' for i = 0, 1, . . . , 2^ 

i = 0,l,...,2'=-l} 

where k < k{n) . 

P{Gl U HI) < 2^(2'' + l)X~^f^2-^'' + 2'=(2'= + l)A-'=^2-'='^ 



by Markov inequality and (4.1) as k < k(n) 

< c(2(2--)a-^)'= 

= C7*= where 7 = 2^^-''^-'^ < 1. 
Given e > and rj > 0, choose ko such that 

c ^ 7*'' < ?7 and 4 ^ A^ < e. 

k>ko k>ko 

Choose 5 = 2~^° . If 5 < 5„ for some n > no, then wx„{^) < wx„{Sn) and we have 
P(wXn{S) > e) < r/ by (4.2). A little more work is required to show that (4.3) holds in the 
< <J) i-e. fco < k{n), case. Pick any (s, q), (t, r) G x (^nN where 5n = 2~^^"^\ such that 
|s— 1|, \q—r\ < 6. We can find a sequence of points (s, q) = (si, qi), (s2, q2)---{sm,qm) = {t, r) 
(refer to lemma 4.1) such that on the event nfco<fc<fc(n)(^'= ^ ^k) we have 

m—l 

\Xn{s,q) - Xn{t,r)\ < \Xn{si,qi) - Xn{s,+i,qi+i)\ < 4 ^ A^ < e. 

1=1 k>ko 
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Now 



Therefore, 



ko<k<k{n) ko<k<k{n) 

Piw^JS) < e) > P{ fl {GknHk))>l-rj. 

ko<k<k{n) 



Thus (4.3) is satisfied with 6 = 2-^o. 



□ 



We now apply Proposition 4.1 to processes X„ = n ^/^Yn{t, r) where Yn{t, r) 
EYn{t,r). 

4.1 Verifying the first tightness condition 



Ynit,r)- 



We check that (4.1) holds for n ^^'^Yn{t,r). Let a > and 

5/4 + a<p < 3/2. 



(4.4) 



We show that there exist constants a > 2 and < (7 < oo independent of n, such that with 
6n = n~^, for all n (sufficiently large) 



E (n-'/'\Yn{s,q)-Ynit,r)\y^ < C\{s,q) - (t,r)r 



(4.5) 



for all is,q),{t,r) G [0,T] x [-S,S] with \{s,q) - (t,r)| > Sn- 

We can assume WLOG that s < t in the following calculations. Define 



and 



Am,j = {Xm,j{ns) < [q^/n ] + [nvs], Xmjint) > [r^/n ] + [nvt]] 
Bm,j = {Xm,j{ns) > [q^/n ] + [nvs], Xm,j{nt) < [r^/n ] + [nvt]} . 
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If g < r then 



r)J{m) 

Yn{s, q) - Yn{t, r)= ^ l{Xm,j{ns) < [qV^ + [nvs]} 

+ ^ ^ l{^mj(ns) < [q^/n\ + [n7;s]} 

m=[ijv^]+l j=l 

- ^ ^ l{Xmj(ns) > [qVn] + [nvs]} 

m<[q^] i = l 

- Yl ^i^rnjint) < [r^/n\ + [nvt]} 

m>[rVn] j=l 

+ l{^m,,(nt) > [r^] + [nvt]} 

m=[gVn]+l j=l 
7)5' (m) 



(4.6) 
(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 



Combining (4.6) and (4.9), (4.7) and (4.10), (4.8) and (4.11) and adding and subtracting 

l{Xmj(ns) > [gV^] + [nvs],Xm,j{nt) < [r^/n] + [nvt]} 

m=[(jv^]+l i=l 



we get 



where 



Yn{s,q)-Yn{t,r) 



[rVn ] 

meZ m=[(jv^ ]+l 



Gm= Y (1^™. - IS'-J - PoMiPiAn.,l) - P{Bm,l)). 
j=0 



(4.12) 



Similarly, when q > r 

Yn{s,q) -Ya{t,r) 



[qy/n ] 
m=[ry'n ]+l 
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Using the identity (a + 6)'^ < 2^{a^ + we get 



E (n-i/^(y„(s, q) - Yn{t, r))) < ^ {EA^' + EB^^) (4.13) 



where 



and 



m 



[(rV(7)v^ ] 
m=[{r/\q)\/n ]+l 

We now bound ElA^"^] and ^[^^^J. 

To bound EB^^, we use the following lemma which is a slight modification of Lemma 8 
in do]. 

LEMMA 4.2. Let 1^ be independent random variables with E[\Yi\'^^] < c < oo and EYi = 

for all i and for some fixed r > 0. There is a constant C < oo such that, for any n , 



^1 (n + >2 + • • • + < C{2r)\rf. 



Proof. Since EYi = 0, 

e[{y. + y, + ... + y„) = E V^,,^'^',j ^^"^^" • • • E-^n- 

where extends over all n-tuples of integers ri, r2, . . . r„ > such that each / 1 and 
ri H h r„ = 2r. 

|£'Y']^^£^y2^^ ■ ■ ■ -^^"1 ^ c by the bounded moment assumption and Holder's inequality. 
The number A = ^^'(ri! • • • r^!)"^ is the coefficient of x^"^ in 

and consequently A < "^"If ^ for x > 0. But y'7>o ^ < 1 + j;^ for x < 1, so taking x = n^^/^ 

gives ^ < + < C'"''- □ 

Recall that viq{x),x G Z are independent with mean Pq{x). Let C denote a constant 
that varies from line to line in the string of inequalities below. Applying Lemma 4.2 to the 
term in equation (4.13) with r = 6 and by the moment assumption on r/Q, we get 

^[5^2] < C(12)!|[r^] - [qVn\f < C{\r - g|V + 1). (4.14) 

We use the following lemma to bound £^[yl"'^^]. 
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LEMMA 4.3. There exists a constant C such that for each positive integer 1 < k < 12 
and for all m, 

E\\Gm\ ] < Dm 

where Dm = C{P{Am,i) + P{Bm,i)]. 

Proof. The proof is the same as in Lemma 4.7 (pages 784-785) of [13]. □ 
LEMMA 4.4. There exists a c > such that 

E[A^^]<c\^+[Y.Dryj |. (4.15) 

Proof. A = Ylmei^rn- Note that EGm = for ah m and the GmS are independent. Let 
X]^'^^ denote the sum over ah /c-tuples of integers ri, r2, . . . , > 2 such that r\ + ■ ■ ■ + rk = 
12. 

mi,...,mi2eZ 

^EE^'^,:^^ E E\GmrE\Gmr---E\GmJ^ (4.16) 

k=l m\^m2^---mj^ 

^ E E E sioj- ■ ■ ■ E '^IG-.I'-' 

k=l i- ^- m m m 

Since < Dm for all m G Z and 1 < / < 12 we get, for aU 1 < A; < 6, 

E\Gmr E ^1^™!''' • • • E ^l^'mf^''^ < max{l, ( Z)„)6} < {1 + (^ D™)^}. 

m m m m m 

Thus 

<c|l+ ^J^D™^ |. (4.17) 

□ 

We evaluate X^.^^^ Dm below. Recall that 



meZ VmeZ meZ / 
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^ P{Am,i) = ^ P{X{ns) < [qy/n] + [nvs] - m,X{nt) > [ry/n] + [nvt] 
^{P(X(ns) = [qV^] + [nvs] - I) 



m) 



m&Z 



meZ l>r 

X P[X{nt) — X{ns) > [nvt] — [nvs] + [r\/n] — [(7\/n] + ^ ~ ''^)} 
= E T^^PiXins) = [qV^] + [nvs] - I) 

l& k>0 

X P{X{n{t - s)) > [nvt] - [nvs] + [r^/n] - [g\/n] + k)} 

< ^ P{X{n{t - s)) - [nv{t - s)] > [r^/n] - [q^/n] + k) 
k>0 



(4.18) 



Similarly, 



P{Bm,i) < P{X{n{t - s)) - [nv{t - s)] < [r^^ - [qV^ + k + 1) (4.19) 

m£Z k<0 



Together, 



P{Am,i) + Y PiPrn,i) < E\X{n{t - s)) - [nv{t - s)] + [qV^] - [rV^]\ 

< E\X{n{t - s)) - [nv{t - s)] \ + \r - q\y/n + 1 

< c{^/n{t — s) + \r — q\\/n + 1}. 



Consequently, (4.15) becomes 

E[A^^] < c{l + {n{t - s)f + \r- q\^n^}. 



(4.20) 



Putting ( |4.14[ ) and ( |4.20[ ) together in ( |4.13| ) we get 



E[{n-'/^{%{s, q) - y„(t, r))}i2] < c{n-' + (^t - sf + \r - qf) 

Using P < 3/2, |r — g| < 2S < oo and t — s < T < oo, we can find constants c > and 
a > 2 such that 

E[{n-'/HYn{s, q) - Yn{t, r))}'^ ] < c{ \r - q^ + \t - sD, 



if \r — q\ > n ^ oi t — s > n ^ . This verifies the first tightness condition (4.1 ). 
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4.2 Verifying the second tightness condition 



To verify (4.2) for n ^^'^Yn(t,r), it is sufficient to show 
LEMMA 4.5. For any < T, S < oo and e > 0, 

\Yn{t,r) - y„(n-^fci,n~%)| > n'/^ 



lim P 

n— >oo 



u 



0<ki<[Tnl^] 



sup 

fcl7i-''<t<(fci+l)n-'3 



[-n'^S]<k2<lnf^S] '=2'" ''<r-<(fc2+l)n- 



(4.21) 



Proof. Recall from (4.4) that 5/4 + a < P < 3/2, a > 0. We first show that particles 
starting at a distance of n^/'^+°' or more from the interval ([— (5* + l)-y/n], [S-^/n]) do not 
contribute to l^(-,r) during time interval [0,T], r G [—S, S], in the n — > oo limit. 



LEMMA 4.6. Let 

Ni = 



(4.22) 



J2 Yl ^{^nrjint) >-[{S + l)V^ + [nvt] 

m<h(S'+l)v^-nV2+ti j=l 

for some < t < T} 
for some < t <T}. 

Then ENi ^0 as n ^ oo. 

Proof. Choose a positive integer M large enough so that 1/2 — a{2M — 1) < 0. The 
expectation of the first sum in (4.22) is bounded by 

C 



E 



P{ sup {X{nt) - nvt) > -[{S + 1) V^] - 1 - m) 

0<t<T 



< C y P( sup {X{nt) - nvt) > I - 1) 
<C l'^^'E[{X{nT)-nvT)l^] 

;>„l/2+a 

by application of Doob's inequality to the martingale X{t) — vt 



<c Y 

/>ni/2+a 
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as E[{X{nT) - nvT)f^] is O(n^0 

< Cni/2-"(2A^-i) ^ as n ^ oo. 
Similarly the expectation of the other sum goes to as n ^ oo. □ 
Let 

([5yri])+nl/2+a. 

m=[-(S+l)v^]-nl/2+"; 
meZ 

be the number of particles initially within n^f^+o^ distance of the interval {—[{S+l)^/n], [S^/n]). 
Fix a constant c so that 

lim P{N2 > cn^^^+") = 0. (4.23) 

n— »oo 

Consider the event 

(j I sup |y„(t,r) -y„(n-%,n-%)| >ni/4e|. (4.24) 

0<fci<[rn''] Uin-'3<t<(fci+l)n-'' J 
[-n'35]<fc2<[n'35] k2n-^<r<{k2+l)n-l' 

If = n~^ki and rg = n~^k2, then 

|^n(t,r) -y„(to,ro)| < |y„(t,r) -y„(to,r)| + |y„(io,0 -^n(to,ro)|. 



For fixed ki and A:2, the event in braces in (4.24) is contained in the following union of two 
events: 



sup sup 

to<t<to+n~f^ ro<r<ro+n 



Ynit,r)-Yn{to,r)\ > \en^'^^ (4.25a) 
U| sup |yn(to,r)-y„(to,ro)| > ^enVn (4.25b) 



The first event (4.25a I implies that at least one of the following two things happen: 



1. At least ^ev}/^ particles cross the discretized characteristic 

s ^ [r\/n] + [nvs] 



for some r E [ro, ro + n during time interval s E [to, {to+n ^)] by jumping. On the 
event {A'"i = 0}, these particles must be among the N2 particles initially within 
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Fi gure 1: Two characteristic lines at distance 5-^ 



fn = n ^ \fn apart. 



distance of the interval (— [(5* + 1)\/?t]; Therefore, conditioned on {A''! = 0}, 

the probability of this event is bounded by the probability that independent rate 
1 random walks altogether experience at least ^en^/'^ jumps in time interval of length 



n 



1-/3 



2. At least \en}l'^ particles cross the discretized characteristic 

s ^ [r\/nl + \nvs\ 



for some r G \rQ^rQ + n~^], during time interval [io, (^o + by staying put while 

the characteristic crosses the location of these particles. These particles must lie in 
the interval J at time nto, where 



J 



[roV^ ] + Mo], [(ro + n ] + {nv{tQ + n 



(refer to figure [T]). For large enough n, the distance between the endpoints of J is at 
most 2. So these \en}l'^ particles must sit on at most 2 sites, say and 



The second event (4.25b) implies that at least \en}-l^ particles either lie in the interval 
[[rov^] + ["'^'io]) [r-y/rij + [refto]] at time nto, or lie in the interval [[^'ov^j ['^x/"]] at time 
0. Since 
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and 

[[ro\/n], [r^n]] C / 

where / = [[ro-^/n ], [{tq + n~^)^/n ]] (refer to figure tliis event implies that at least |en^/^ 
particles lie in interval J at time nto or at least ^en^/^ particles lie in interval I at time 
0. For large enough n, the distance between the endpoints of I is at most 1, so the |en^/^ 
particles lying in interval / at time must sit on a unique site Xq say. 

Let n(cn'^/^^"~^) denote a mean cn^/'^~^°'~^ Poisson random variable that represents 
the total number of jumps among cn^/^^° independent particles during a time interval of 
length n^~^. Then, 



P( event in ( |4.2lD ) 
< P (iVi > 1) + P (^N2 > ' 

[Tn/3] [SnP] 

+ E \P (n(cn3/2+-/3) > l^„l/4^) + p ^^.(^0 > l^^l/4^| (4.26) 



A;i=0 k2 = l-Snf^ 



4.6 



and 



The probabihties P{Ni > 1) and P{N2 > cn^/^+a) vanish as n ^ oo by lemma 
(4.23). n(cn^/^+"~^) is stochastically bounded by a sum of Mn = [m3/2+«-/3] + 1 i.i.d. 
mean 1 Poisson variables, and so a standard large deviation estimate gives 

P (n(m3/2+"-/3) > len^/^) < exp {-MnI{lM-'en'/^)] , 

where / is the Cramer rate function for the Poisson(l) distribution. By the choice of a and 
/5, Mn > n", while M~^n^/^ oo. Consequently, there are constants < Co, Ci < cxd such 
that 

^ ^ P (n(cn3/2+"-/3) > ien^/^) < Con^^ exp{-Cin-} ^ 0. 

A:i=0 fc2=[-Sn'3] 

By Lemma 4.10 in [13], we have 



So, 



sup £;[?7i(x)^^] < oo. (4.27) 

x(^Z,t>0 



fci=0 fe2 = [-5n/'] 

< {Tn^ + l)(25n^ + 1)8^2^-12^-3 E['n'^{x)^^] 

x,t,n 



29 



vanishes as n ^ oo by (4.27) and because 2/? — 3 < 
Similarly for the other probability in (4.26). 



□ 



Since the two conditions of Theorem 
we can conclude that 



4.1 



hold for the sequence of processes {n ^^^Yn}, 

(4.28) 



limlimsupP{tt;„-i/4y (S) > e} = for all e > 0. 

|5J,0 n " 



4.3 Weak Convergence 

Finally, we use the theorem about weak convergence in D2 from [2,. By Theorem 2 in [2] 
we have X„ converges weakly to X in D2 if, 



1. (X„(ti,ri),--- ,Xn{tN,rN)) converges weakly to {X{ti,ri), 
nite subsets {{ti,ri)} £ [0,T] x [— S, S], and 

2. lim^^o Ihiisup^ > e} = for all e > 0, where 

Wx{5) = sup \x{s,q) - x{t,r)\. 

(s,q),{t,r)(^[0,T]x[~S,S] 
\{s,q)~{t,r)\<5 



I X{t]\f, rj^f)) for all fi- 



This, together with the convergence of finite-dimensional distributions of {n~^^^Yn{-, •)} 
and (4.28) gives us weak convergence of {n''^^^Yn(-, •)} as n — > 00. Since the expectations 
n~^^^EYn{t, s) vanish uniformly over < t < T , < \s\ < S hy (3.4), we conclude that the 
process {n~^^^Yn{-, •)} converges weakly as n — > 00. 



5 Proof of large deviation results 

Proof of Theorem 2.2 and Corollary |2.2[ Assume that r]Q{m),m G Z are i.i.d. Fix 
r G M and t > O.We prove that n^^^'^Yn{t,r) satisfies the LDP with a good rate function. 
We start with some preliminary calculations. 

Yn{t,r)= 5^ [l{^m,i(nt) < M + [rv/^]}l{m > [rV^]} 

m=—OD j=l 

- l{Xm,j{nt) > [nvt] + [rVn]}l{m < [r^/r^]}] (5.1) 

m=— 00 j=l 
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Define 



(5.2) 



If m > [r-y/ri] then 



l + (e^-l)i?[/;;«(t,r)] 



(A/:J^(t,r))' 



A:! 



(5.3) 



1 + (e^ - l)P(X(nt) < [nvt] + [r^/n] - m 



where X(-) represents a random walk with rates p{x) starting at the origin. Similarly, if 
m < [r^/n] then 



= 1 + (e"-*^ - l)P{X{nt) > [nvt] + [r^/n\ - m). 
We now calculate the logarithmic moment generating function for l^(t, r). 



fc>0 



(-A-/:f (t,r)) 



kl 



(5.4) 



|m-[rv^|<[nt(5] [ i=l 

i=i 



(By large deviation bounds on X{nt) (3.1) and (3.2) , the second term is of o{^/n)) 



logi?exp<' A ^ [Ci'^(t,r 



f:::fit,r)]}+o{v^) 



\m—[ry/n]\<[nt5] 
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lim ^ log ^e^'^" = lim V logE 

J — ^nn , /r) n — >00 a/ 7X ^ — ^ 



|m— ry^] |< [nt5] 



AE:!r'[c,',^'(*,o-c,?'(*,^)]' 



e 



lim 



\m-r^]\<[ntS] k>0 

(Recall that 7(a) = log Ee'^'^o(-)) 

= lim ^ 7(logM-(A)). 



|m— ry'n] |< 



(5.5) 



Using (5.3), (5.4), (5.5), applying Central Limit Theorem and a Riemann sum argument 

we get. 



A(A) = lim ^ log £;e^^" 



n— >oo 

(•00 



poo pv 

/ 7(log{l + (e^ - l)^,,t{-x)})dx + / 7(log{l + (e-^ - l)<^^,t{x)})dx. 

Jo J-00 

(5.6) 



By Assumption 2.3 we get A(A) < 00 for A G R. It is also easy to check that A(A) 
is strictly convex and essentially smooth on R. By Theorem 2.3.6 in [5] (Gartner-Ellis 
theorem) /(•), the convex dual of A(A), is the good rate function. We now find the explicit 
expression for the rate function. 

If rjQ^m) ~ Poisson(/?), then 7(a) = /o(e" — 1). Therefore, 

poo pO 

A(A) = p{e^ - 1) / <^^^ti-x)dx + p{e-^ - 1) / ^^^t{x)dx 



The convex dual of this is: 
I{x) = sup{xA - A(A)} 



xlog 



+ 1 1 + 



/ x^/tt 



^r. 1 



vr 



1 I for X € 
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This proves (2.15). 



To prove Theorem 2.2 we first check that /(•) given by (2.14) is the convex dual of A(-) 



and then outline how to get the expression in (2.14). Elementary but tedious computations 



give 



/'(A'(A)) = A. 



(5.7) 



It can be shown that A' is continuous and strictly increasing. Therefore, /' is defined on the 



whole real line by (5.7). By Theorem 26.5 in [12] (page 258), I must be A* plus a constant. 
But /(O) = = A*(0), so / = A*. This proves that I given by (|2.14[) is the convex dual 



of A and hence the rate function. Since A' is continuous and strictly increasing, /' must 



be strictly increasing from (5.7) and hence / must be strictly convex. This completes the 



proof of Theorem 2.2 Corollary 2.2 comes as a special case of (2.14). 



For the reader's benefit, here is an indication of how the expression in (2.14) is derived 



non-rigorously. We first approximate the integral in A(-) by a Riemann sum. 



where 



AiiX) :-- 



(A)) 



\57(Sr$^^^(fc5)(-A)) 



for A; > 
for A; < 0. 



(5.^ 



The function Brp{\) denotes the log moment generating function for Bernoulli random 



variables as defined in section 2.3. Observe that the summands (5.8) are a composition of 



two functions. Using the definition of convex dual, it is easy to prove the identity 

(/ o gTix) = ng{\))g*{g'{X)) + nfigm 



(5. 



where A is such that 



x = ng{X))-g'iX). 



(5.10) 



We use this identity to get the convex dual of the summands (5.8) 
For small 6 

A(A)« ^i(^)- 

\k\<ll/S] 



33 



The convex dual of the sum is then given as an infimal convolution. 

A*{x) R 



Jfc|<[l/'5] 



{x)= inf _ J2 (^i) (^'^) 



now using (5.9) we get 



r iv<5] 

I k=l 



inf 

|fe|<[l/5] 



+ 7* (7'(^n-$«2t(fc5)(Afc))^ 



fc=-[l/5] 

+7* (y (^''*.2t(fc5)(-Afc)) 



where 



Xk 



5i (-Bri_$^^^(fc5)(Afc)) (fc5)(Afc) for > 



H [Br^.A'^S){-^k)) Br',J-Xk) for k < 0. 



Note that 



Bri-^^,^,(y){X) = Zx{y) for y > , 5r$^^^(j^)(-A) = Zx{y) for y < 
^^'i-<l..,,fe)(A) = 1 - Fx{y) for y > , ^(^)(-A) = Fy^{y) for y < 0. 



(5.11) 



(5.12) 
(5.13) 



Taking the limit as (5 ^ of (5.11) we get the constrained variational problem 

l-OO 

A*(x)= inf / 7*{y(Z;,(,)(y))}dy 



i^- am b{Zx{y){y))]dy = x} 



oo 
oo 



+ 



7'(^A(y) (y) ) (FA(y) (y) ) dy 



Solving the above variational problem using standard functional analysis techniques we get 
A(y) = a{x) minimizes the above functional and a(x) satisfies the constraint (2.11). □ 
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Proof of Theorem 12.31 Let Yn{t) = Yn{t, 0). Under the assumption rjoim) ~ Poisson{p) 
we can show that {n~^^'^Yn{-)^ satisfies the large deviation principle in /^^[O,^). 
Fix k time points < ti < t2 < ■ ■ ■ < ifc- Define the fc-vectors with 0, 1 entries by 

Fm,j = {l{Xmj{nti) < [nvti]}, l{Xm,j{ntk) < [nvtk]}) 

and 

Gm,j = {l{Xm,j{nti) > [nvti]}, l{Xm,j{ntk) > [nvtk]}). 

Then 

oo ^S{m) Voim) 

{Yn{h),...,Yn{tk)) = J2 E ^-J- E E 
m=l j=l m=—oo j=l 

Let u := {u^^\ . . . , u^''^) E {0, 1}^ be a fc-vector with u ^ 6. Define 

^n(^):=E E HFm,,=u} 

771=1 j = l 

and 

Voim) 
m=— oo j=l 

If 7?o(') i-i-d. Poisson(p) random variables, then N^{u) is a Poisson random variable 
with rate 

oo k 

5^^p(n<-)<oo, 

7Tl=l i=l 

where ^ = {X{nti) < [nvti] — m}, ^ = {X(ntj) > [nvti] — m}, and N'^{u) is a Possion 
random variable with rate 

k 

Y pp(rK^)<oo, 

m=— oo i=l 

where D^- = {X{nti) > [nvti] — m}, D^- = {X{nti) < [nvti] — m}. We use the bounds 
(3.1) and (3.2) on the large deviations of random walks to justify the rates being finite. 
We can write 

. . .,Yn{tk)) = YjN^iu) - K{u))u. 

Let Uj = {u"p , . . . ,u"l'''),j = 1,...,2'^ — 1 denote the {0,l}-valued fe-vectors, excluding 
the zero vector. By the contraction principle (Theorem 4.2.1 in [5]) we can conclude that 
(Yniti), . . . ,Yn{tk)) satisfies the LDP with good rate function given by 

lH,...,tk{x) := inf J(yi, . . .,y2k_]_,zi, . . .,Z2k_i) 
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for any x = (xi, . . . , Xk). The inf here is taken over the set {(yi, . . . , ^2^-1; ^i, . . . , Z2k_i) : 
X = Ylu^o^yj ~ -^jO^j}- J{yi-> ■ ■ ■ ) 2/2^-1) -21, . . • ,-22*-!) is the good rate function for the 
sequence of vectors {{N^{ui), . . . , A^^(?22*_i), iV^(nl), . . . , N!^{u2k_i))}n of 2^~^^ - 2 inde- 
pendent Poisson random variables. Define for x >0, 

= {B{K2t,) < -X}, = {B{K2t,) > -X}, 

DIj = {B{K2tj) > x} and D'^ j = {B{K2tj) < x}, 
where B{-) is standard Brownian motion. 

2'' -I 

J{yi, ■ ■ ■,y2k-i,zi, . . .,Z2k_i) = V [yjlog — + Zilog^ - ai{— - 1) - Pi{^ - 1)] 

^ ai Pi Oi Pi 



where 



i=l 



and 



A = lim -^EN^iiii) = p 

n^oo a/ 77, 



We now apply Theorem 4.30 in to {n ^^'^Yn{-)}. This gives us the large deviation 
principle for {n~^^'^Yn{-)} in /^^[OjOo) with good rate function 

I{x) = sup/ii,...,i„(x(ti), . . .,x{tm)). (5.14) 

□ 
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